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ABSTRACT 

The classical dN/dz test allows the determination of fundamental cosmological parameters from the 
evolution of the cosmic volume element. This test is applied by measuring the redshift distribution of 
a tracer whose evolution in number density is known. In the past, ordinary galaxies have been used as 
such a tracer; however, in the absence of a complete theory of galaxy formation, that method is fraught 
with difficulties. In this paper, we propose studying instead the evolution of the apparent numbers of 
dark matter halos as a function of their circular velocity, observable via the linewidths or rotation speeds 
of visible galaxies. Upcoming redshift surveys will allow the linewidth distribution of galaxies to be 
determined at both z ~ 1 and the present day. In the course of studying this test, we have devised 
a rapid, improved semi-analytic method for calculating the circular velocity distribution of dark halos 
based upon the analytic mass function of Sheth et al. (1999) and the formation time distribution of Lacey 
& Cole (1993). We find that if selection effects are well-controlled and minimal external constraints are 
applied, the planned DEEP Redshift Survey could allow the measurement of the cosmic equation-of-state 
parameter w to ±10% (as little as 3% if fi m has been well-determined from other observations). This 
type of test has the potential also to provide a constraint on any evolution of w such as that predicted 
by "tracker" models. 



1. INTRODUCTION 

For almost as long as they have been known to be lo- 
cated outside the Milky Way, galax ies have been used to 
probe the geometry of the universe ( Hubble 1926| ) . Fairly 
early in the history of relativistic cosmology, it was real- 
ized that if the number density of galaxies is known, their 
apparent numbers per solid angle and redshift int erval de- 
pend directly on basic cosmological parameters ( Tolman] 
1934 ). Applying this u dN/dz test" requires a substantial 
sample of distant galaxies with measured redshifts; it was 
first applied by Loh and Spillar (1986) using a set of 406 
galaxies with estimated photometric redshifts z p between 
0.15 and 0.85. Assuming that the comoving number den- 
sity of galaxies is constant and that their luminosity func- 
tion retains similar shape over that redshift range, they 
measured Q, = 0.9±g*J. 

However, it remains unclear whether the assumptions 
of Loh and Spillar hold; no complete theory or si mula- 



tion of galaxy fo rmation and evolution as yet exists (Ben- 
son et al. 1999). It would thus be preferable to apply 
such a test to objects whose abundance can be calculated 
semi-analytically or via computer models, eliminating such 
ambiguities. In this paper, we suggest that the observed 
numbers of galaxies (and thus indirectly dark halos) as 
a function of linewidth and redshift provide a candidate 
for improving on the analysis of Loh & Spillar, particu- 
larly because such quantities will be measured by surveys 
occurring in the next decade. 

In the process of studying this possibility, we have de- 
veloped a semi-analytic method for determining the abun- 
dance of dark halos as a function of their circular velocity 
= GM/r) and redshift that can provide results similar 
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to those of N-body simulations in considerably less com- 
puter time and without the effects of limited resolution. 
In essence, we begin with the abundances of galaxies as a 
function of mass predicted by the modified Press-Schechter 
formalism of Sheth et al. (1999), which provides an excel- 
lent fit to dark-matter simulations. We then determine 
the probability that a halo of given mass has formation 
time such that it will virialize with the desired circular ve- 
locity using the semi-analytic method of Lacey and Cole 
(1994), which again is very well reproduced by simulations. 
By integrating over the distribution in mass we may then 
determine the total abundance of objects with a given cir- 
cular velocity at the epoch of interest. We describe this 
method in more detail in §2, and the prospects for our 
proposed test in §3. 

2. DERIVING THE VELOCITY FUNCTION OF DARK HALOS 

A number of semi-analytic methods for calculating the 
veloctity d istribution of dark halos based upon the Press- 
Schechter (Press & Schechter 1974) formalism have been 
proposed in the past. Narayan & White (1988) used the 
assumption that halos have just virialized at the epoch of 
observation to map the Press-Schechter distribution of ha- 
los in mass into a distribution in velocity for a f2 = 1 Cold 
Dark Matter (CDM) universe; Kochanek (1995) extended 
their analysis to other common cosmological models. How- 
ever, the assumption that objects observed are only just 
virializing is clearly inappropriate in many cases, particu- 
larly for lower-mass halos; once objects virialize and break 
away from the cosmological expansion, their circular ve- 
locities evolve little if at all save in major mergers (Ghigna 
et al. 1999). To account for this, Kitayama and Suto 
(1996a) proposed that the distribution of X-ray temper- 
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atures of clusters (T oc t£) be calculated by totalling the 
rate at which objects of the desired temperature that sur- 
vive without merging into a much more massive object 
through the epoch of interest form, using the results of 
Lacey & Cole (1993). As the authors of that paper noted, 
there remain some ambiguities in their formulation, how- 
ever. 

We have attempted to improve upon these prior meth- 
ods for calculating the velocity distribution of dark matter 
halos, taking advantage of the progress made over the past 
few years in extending the Press-Schechter formalism. It 
has been repeatedly demonstrated that standard Press- 
Schechter calculations overpredict the abundance of low- 
mass halos compared to equivalent simulations; a number 
of remedies have been suggested recently (e.g. Bond & 
Myers 1996). One of the simplest to implement is that of 
Sheth et al. (1999), which accounts for the average ellip- 
ticities of dark matter halos by making the critical over- 
density for their collapse mass-dependent. The resulting 
analytic expression for the number density of halos of given 
mass (given by their equations 1 and 6) provides an ex- 
cellent fit to the results of a variety of cosmological sim- 
ulations. The results of Sheth et al. allow a prediction of 
halo abundances much improved over the standard Press- 
Schechter formula with only minimal complications to the 
calculation. We have therefore used their equations to de- 
termine the distribution of halos in mass at a given rcdshift 
in our work. 

Because the density of the universe was greater at early 
times, objects that virialized at high redshift are more 
compact than those that virialized more recently, and 
therefore have a higher circular velocity for their mass. 
Furthermore, once it has collapsed, the circular velocity of 
a halo changes only minimally (Ghigna et al. 1999). There- 
fore, in determining the abundance of galaxy-scale halos 
with a given circular velocity today, it is preferable to take 
into account the time of formation of those halos in addi- 
tion to their mass. Fortunately, this may be calculated; 
Lacey & Cole (1993) derived a semi-analytic formula for 
the distribution of formation times of halos with a given 
mass that exist as virialized objects at a given redshift 
(defining the formation time as the epoch when the object 



first h |ad at least half its final mass). This formula provides 
an excellent tit to the results ot simulations (Lacey & Cole 
1994), even though it assumes spherical collapse while the 
simulated halos could be ellipsoidal. We therefore have 
used their semi-analytic formula (equation 2.19 of Lacey 
& Cole 1994) to determine dp/dif, the probability dis- 
tribution of the epoch of collapse if of halos observed to 
have mass M at time t a . Note that the integral defining 
dp/dtf must be evaluated with particular care where the 
integration variable M\ sa M . However, an approximate 
analytic solution exists in that region, and the integral 
may be safely performed over the remainder of its range 
numerically. 

In practice, we are interested in calculating the abun- 
dance of halos at redshift z that have circular velocity 
v c ; for each halo mass, we therefore need only evaluate 
dp/dtf at that formation epoch for which a halo of the 
given mass would virialize with the desired circular ve- 
locity. We may obtain that time from the relationship 
between the circular velocity and the mass M, which is 
most conveniently parameterized by using the length scale 



ro (/i _1 Mpc) = (3M/47rpo) 1 ^ 3 , where po is the present-day 
mass density of the universe. By combining the definitions 
of circular velocity, of ro, of A V i r (the ratio of the mean 
density of a virialized halo to the critical density of the 
universe when it virialized), and of the dimcnsionlcss evo- 
lution of the Hubble parameter E{z) we find (cf. Narayan 
& White 1988, Kochanek 1995): 
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where we have determined the circular velocity for a halo 
of mass M/2, more appropriate given our definition of the 
formation time. The objection might be raised that we 
are using a formula for circular velocity that implicitly 
assumes spherical geometry but are applying it to poten- 
tially ellipsoidal dark matter halos; however, Lacey & Cole 
(1996) found that the velocity dispersions of dark matter 
halos calculated using a singular isothermal sphere model 
(for which a = \/2vc) matched the properties of N-body 
simulated halos quite well. In our calculations we have 
used the fitting formulae for A V i r given by Bryan and Nor- 
man (1997). 

Thus, to calculate the comoving abundance of halos with 
a given circular velocity at an epoch of interest, we simply 
calculate the distribution of halos in mass on a 500-element 
grid, determine the probability that those halos have the 
desired circular velocity, and then integrate numerically 
over that grid. This procedure takes ~ 5 sec on a modern 
workstation to calculate the abundance at a single circu- 
lar velocity and epoch, far faster than simulations could be 
run, and has none of the limitations of resolution or sample 
size that affect even the largest N-body simulations today. 
Our technique seems to be similar to one employed by 
Kitayama and Suto (1996b) to check their methods, save 
for the use of the ellipsoidal-collapse mass function rather 
than a traditional Press-Schechter calculation; however, 
details on the latter are minimal. 

Our procedures could certainly be improved upon. For 
instance, the connection between circular velocity, forma- 
tion time, and mass could be made using a more realistic 



model than a singular isothermal sphere (Navarro, Frcnk 
& White 1995|, Bullock et al. 1999). Calculation of forma- 



tion times using an ellipsoidal-halo abundance paradigm 
rather than the Press-Schechter method of Cole and Lacey 
might also yield modest improvements. Finally, we must 
note that the mass function of Sheth et al. was based upon 
a fit to mass functions of halos identified in N-body simu- 
lations using the "friends-of-friends" algorithm with selec- 
tion parameter b — 0.2; in such a method, as in standard 
Press-Schechter calculations, smaller subhaloes within a 
more massive halo (e.g. galaxies within a cluster) are not 
counted separately in the mass function. Unless an im- 
proved method taking this into account is used, we can 
only calculate the abundance of relatively isolated dark 
matter halos (i.e., those that are not part of a larger, viri- 
alized group or cluster). Such objects are readily identified 
observationally in well-sampled redshift surveys, but this 
is still a significant limitation. Despite these caveats, we 
believe that our technique for semi-analytic calculations 
provides a substantial advance over methods previously 
used; all of them suffer from similar flaws and have fewer 
advantages. Our methods may also be useful for a variety 
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of cosmological calculations; for instance, the circular ve- 
locity distribution of dark halos is required for calculations 
of the abundance of strong gravitational lenses, and hence 
for the limits placed on fl^ thereby. 

3. THE PROPOSED TEST 

In order to study the evolution of the volume element 
of the universe, we require some identifiable tracer whose 
number density evolution is well known. The dark halos 
which contain galaxies are an obvious candidate; the re- 
sults of semi-analytic calculati ons (e.g., Shcth et a l. 1999) 
and N-body simulations (e.g., Jenkins et al. 1997) for the 
behavior of dark matter are by now fairly well understood. 
Of course, directly observing the properties of dark ha- 
los at high redshift is impossible; however, some of their 
parameters may be determined by studying galaxies ly- 
ing within them. In particular, it is well-known that the 
motion of material in the outer parts of typical galaxies 
is dominated by the gravitational influence of dark mat- 
ter (Fischer et al. 1999); the depth of the halo potential 
well may thus be directly derived from observed linewidths 
or velocity dispersions of galaxies. Unlike mass measure- 
ments, such observations are relatively insensitive to dis- 
tance from a galaxy's center (for instance, the rotation 
curve s of many spiral galaxies are observed to be fairly flat 



pleted. This survey, to be begun in early 2001 and com- 
pleted within the five years following, will obtain spectra of 
~ 60,000 galaxies at redshift z > 0.7 with a resolution of 
~ 80 kms -1 FWHM, allowing measurements of linewidths 
for perhaps a quarter to half of the galaxies in the sample. 
The number of spectra obtained will be great enough to 
test for selection effects in redshift or linewidth by com- 
parison to calculations or simulations at the few percent 
level. For instance, smaller halos are predicted to be dis- 
tributed in circular velocity according to a power-law with 
a slope well-determined by both N-body simulations and 
our techniques; deviation from such a distribution would 
be a strong indicator of observational selection effects, al- 
lowing us to correct for or avoid regions of parameter space 
suffering from incompleteness. Complementary informa- 
tion on the linewidth distribution of galaxies at low red- 
shift should be availa ble fro m the Sloan Digital Sky Survey 
( |Loveday et al. 1998| ), 2DF flCollcss 98j ), or other local sur- 
veys, allowing accurate normalization to the abundance of 
halos at z = 0. 

The determination of the volume element afforded by 
such measurements can place constraints not only on the 
simplest models of the universe which include only mat- 
ter and a cosmological constant, but also on so-called 



quintessence" models (Turner fc White 1997; Caldwell 



over a | great range of radii, as would be predicted for halos Dave, fc Stcinhardt 1998| ) in which the equation of state 



which are roughly isothermal spheres). Under modest as- 
sumptions (e.g. about their radial profiles), the expected 
abundance of dark matter halos as a function of observable 
linewidth may be determined from simulations or from the 
semi-analytic method described in §2. 

Results of our semi-analytic calculations for the evolu- 
tion of the differential comoving abundance of isolated ha- 
los with circular velocity 200 kms -1 in some commonly- 
used cosmological models are depicted in Fig. 1. A CDM- 
like power spectrum with shape parameter Y = 0.25 and 
normalized to match the value of ag derived by Borgani 
et al. (1999) from X-ray observations of galaxy clusters 
was employed in all cases. As it is likely that selection 
effects will be somewhat similar for galaxies at z — and 
z ~ 1, we plot the ratio of the abundance of galaxies at 
a given redshift to that at z = 0. The abundance of ha- 
los at low circular velocities is smaller now than in the 
recent past because of merging of smaller halos into new, 
larger objects and because smaller halos lose their individ- 
ual identities in a Press-Schechter framework when they 
are included in a larger virialized object such as a group 
or cluster. Qualitatively, the predictions of N-body sim- 
ulations are similar (see, e.g., the cluster simulations of 
Ghigna et al. 1999). Our calculations demonstrate that 
the comoving abundance of halos with fixed linewidth at 
z ~ 1 relative to today is very insensitive to the back- 
ground cosmological model under reasonable assumptions 
about the matter power spectrum; e.g. for models with 
a low matter density £l m = 0.3, the relative abundance 
is only 2% different in a flat model from an open one. 
Thus, observations of the apparent numbers of such ob- 
jects per unit redshift per steradian in comparison to their 
abundance today directly yields a measurement of the cos- 
mological volume element, and hence of fundamental cos- 
mological parameters. 

A measurement of dN/dz dVL will be quite feas ible once 
the DEEP Redshift Survey (|Davis fe Faber 1998j) is com- 



of a non-material component, P = wp, can take arbitrary 
value between -1 and 1 (w = for matter, while w = — 1 
corresponds to a cosmological constant). The amount of 
comoving volume per unit redshift per steradian at z ~ 1 
depends strongly on w, as illustrated in Fig. 2. The test 
we propose in combination with a determination of fl m 
(e.g. from velocity statistics of galaxies) should allow di- 
rect determination of the equation of state of any dark 
energy present, as illustrated in Figure 3. In compari- 
son, CMB measurements planned for the near future and 
ground-based observations of the apparent magnitudes of 
Type la supernovae can determine wto~ 10% at best (Hu 
et al. 1999); similar accuracy might be achieved using the 
stati stics of lensed sources in the Sloan Digital Sky Sur- 
vey ( |Cooray fe Huterer 1999| ). Our estimated constraints 
are similar in strength to those predicted for the proposed 
SNAP satellite (Perlmutter et al. 1999), though less sensi- 
tive to the present-day curvature and more sensitive in the 
orthogonal direction (i.e., more complementary to CMB 
measurements). The technique we propose is at least as 
effective as these even without including complementary 
information from other tests. As Fig. 3 illustrates, since 
fl m will be increasingly constrained from other experi- 
ments in the near future (including other analyses of the 
DEEP dataset), the test we propose has the potential to 
determine w to ~ 3%, even without other complementary 
measurements. If the DEEP survey suffers from incom- 
pleteness, there must be at least as many halos of a given 
linewidth as we observe; thus, a lower limit on the volume 
element, and hence on —w, would be provided. 

The test we propose could also detect evolution of w, 
if any exists and other cosmological parameters are well- 
determined. It has recently been proposed that in rea- 
sonable quintessence models w may vary with the current 
matter density as the universe evolves (so-called "tracker 
solutions" ; cf. Steinhardt et al. 1998) or even oscillate over 
time; such changes would be quite difficult to measure via 
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methods that study integrated quantities (e.g., the lumi- 
nosity distances of SNe la). However, the test proposed 
here allows measurement of a more localized quantity, the 
volume element; hence, it may be able to measure changes 
in that quantity even on relatively small scales. By study- 
ing the apparent abundances of the dark halos of galaxies, 
we may test the most recent cosmological models - or even 
those not yet envisioned - with techniques having their ori- 
gin in the earliest days of extragalactic astronomy. 
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fig. i — The evolution of the differential comoving number density of halos with circular velocity 200 
km s" 1 normalized to the value at z — . We consider three common cosmological models: an open 
Cold Dark Matter (CDM) model with Q = 0.3 today; an Einstein-De Sitter model with Q — 1; and 
a flat LCDM model having Q m = 0.3 today and a cosmological constant. Note the similarity of the 
predictions at z ~ 1. 
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fig. 2.— The comoving volume per unit redshift per steradian, dV/dzdQ, in a variety of cosmological 
models (cf. Peebles 1993). The thickest curves indicate the models considered in Fig. 1; the thick solid 
black curve depicts the volume element in a fl m = 0.3 open universe, the red dot-dashed curve shows an 
Einstein-de Sitter model, and the blue dashed curve depicts an LCDM model. The green, light dotted 
curves show the evolution of the volume element in standard quintessence models with Q m = 0.3, 
Qq = 0.7, and constant w = —0.2, -0.4, -0.6, and -0.8 (listed from lowest- volume to highest). Finally, 
the thin, solid purple line depicts a model in which w is proportional to the expansion parameter, 
w = —0.8a. 
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fig. 3.— Predicted 95% confidence level contours for application of the dN/dz test to the DEEP 
Redshift Survey (solid red curves), along with other existing or predicted constraints. These have 
been obtained by conservatively assuming that the DEEP Redshift Survey will obtain useful linewidth 
information on 10,000 galaxies distributed as (1 + z)~ 2 among 8 bins between z=0.7 and 1.5 (i.e., that 
the great majority of observed galaxies are excluded from the test because of potential systematic 
errors or incompleteness). The top panel depicts confidence contours in the Q m — Q\ plane (i.e., under 
the assumption that w = —1) about the point (Q m = 0.3, fl\ = 0.7). Overplotted are the current 
68% confidence SNela constraint from Perlmutter etal. 1 1999 (black, dashed line), the target 95% 
confidence interval for the proposed SNAP satellite (dot-dashed black curve; taken from figures on 
the SNAPSAT website at |http:/ /snap.lbl.gov| ), and a sample CMB constraint (Melchiorri etal. 1999) 
with errors reduced to simulate a determination of fl CU rvature to ±0.1, representative of what may be 
expected from experiments performed before DEEP is completed such as BOOMERANG and MAP. 
The solid red curves in the bottom panel show the predicted 95% confidence contour from the DEEP 
Redshift Survey in the Q m — w plane (under the assumption of a zero- curvature universe) about the 
points (fl m = 0.3, w = —0.7) and (Q m = 0.3, w = —1). Also plotted are target 95% confidence 
constraints for the latter model from SNAPSAT (black dot-dashed curve) and from observations of 
cluster abundances by the proposed COSMEX satellite (blue dotted region; cf. Haiman et al. 2000). 
The techniques involved (and thus, almost certainly, any systematic errors) differ greatly among these 
methods, making each a valuable check on the others. 



